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ABSTRACT 

We investigate a model of galaxy clusters in which the hot intracluster gas 
is efficiently heated by dynamical friction (DF) of galaxies. We allow for both 
subsonic and supersonic motions of galaxies and use the gravitational drag for- 
mula in a gaseous medium presented by Ostriker (1999). The energy lost by the 
galaxies is either redistributed locally or into a Gaussian centered on the galaxy. 
We find that the condition of hydrostatic equilibrium and strict energy balance 
yields a trivial isothermal solution T iso , independent of radius, or rising temper- 
ature distributions provided T iso /j < T < T iso , where 7 is the adiabatic index of 
the gas. The isothermal temperature corresponds to the usual scaling relation be- 
tween the gas temperatures and the velocity dispersions of galaxies. However the 
minimal temperature associated with the rising solutions is ~ |T v ; r , larger than 
that inferred from observations, the radial distribution of galaxy masses notwith- 
standing. Heating by supersonically moving galaxies cannot suppress thermal 
instability, although it can lengthen the growth time up to the level comparable 
to the ages of clusters when Mach number of galaxies is less than about two. 
We show using numerical hydrodynamic simulations that DF-induced heating is 
generally unable to produce stable equilibrium cores by evolving arbitrary non- 
equilibrium clusters, although it can lengthen the cooling time. We conclude that 
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DF-induced heating alone is an unlikely solution to the cooling flow problem, al- 
though it can still be an important heat supplier, considerably delaying cooling 
catastrophe. We discuss other potential consequences of DF of galaxies in galaxy 
clusters. 

Subject headings: galaxies: clusters, general — galaxies: interactions — galaxies: 
kinematics and dynamics — cooling flows — X-rays: galaxies — instability 

1. Introduction 

Galaxy clusters are the largest gravitationally bound structures in the universe. They 
typically contain several hundreds to thousands of galaxies orbiting in a gravitational po- 
tential well formed primarily by dark matter (e.g., Bahcall 1977). They are also filled with 
hot gas with T ~ 2 — 10 keV that loses thermal energy prolifically by emitting X-rays. In 
the absence of any heat sources, the radiative cooling in the cores of rich clusters would 
result in a cooling flow in which gas settles in the gravitational potential and drops out as 
cold condensations (e.g., Fabian 1994). However, recent high resolution Chandra and XMM- 
Newton observations have not shown the expected signature of gas cooling below one-third 
of the virial temperature (Blanton et al. 2001; Peterson et al. 2001, 2003; Tamura et al. 2001; 
Johnstone et al. 2002; Allen et al. 2001). This strongly suggests that there must be some 
source of heat that manages to balance the radiative cooling and thus prevents the mass 
dropout at the central regions of galaxy clusters. Most theoretical models proposed so far 
for cluster heating fall into the following two groups: (1) energy injection via radiations, jets, 
outflows, bubbles, or sound waves from a central active galactic nucleus (Ciotti & Ostriker 
2001; Churazov et al. 2002; Kaiser & Binney 2003; Ruszkowski et al. 2004; Roychowdhury et 
al. 2004; Vecchia et al. 2004 and references therein); (2) diffusive transport of heat from the 
outer regions of the cluster to the center via conduction (Tucker & Rosner 1983; Bregman 
& David 1988; Narayan & Medvedev 2001; Voigt et al. 2002; Zakamska & Narayan 2003; 
Kim & Narayan 2003a) or turbulent mixing (Cho et al. 2003; Kim & Narayan 2003b; Voigt 
& Fabian 2004; Dennis & Chandran 2005). 

Although less well studied than those mentioned above, there are clearly other energy 
sources in the cluster environments. These include kinetic energy in the orbital motions of 
galaxies (Miller 1986; Just et al. 1990), gravitational potential energy of the gas (Markevitch 
et al. 2001; Fabian 2003), feedback from intracluster supernovae (Domainko et al. 2004), dis- 
sipation of turbulent energy driven by infall of subclusters (Fujita et al. 2004), etc. Amount 
of available energy in each source is comparable to or larger than thermal energy of the 
intracluster medium (ICM), so that radiative cooling would be easily offset if there exist 
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physical mechanisms that can tap the energy from the sources and convert it into thermal 
energy in the ICM. 

One such process is dynamical fiction (DF) which occurs when a galaxy moving around 
a cluster center interacts with its gravitationally induced wake in the ICM (e.g., Dokuchaev 
1964; Hunt 1971; Ruderman & Spiegel 1971; Rephaeli & Salpeter 1980). Because of the 
gravitational drag, the galaxy loses some of its kinetic energy. Turbulence in the ICM 
generated by the superposition of the wakes of many galaxies or by other processes will 
absorb the lost kinetic energy and, in the presence of tangled magnetic fields and/or viscosity, 
turns it into heat at the dissipation scales (e.g., Deiss & Just 1996). Miller (1986) and Just 
et al. (1990) independently estimated the heating rates by DF of galaxies and found that the 
DF-induced heating can fairly well balance radiative cooling of a rich cluster, provided the 
mass-to-light ratio is about 20. Using updated data for the properties of galaxies and gas in 
the Perseus cluster and employing Monte-Carlo approaches, El-Zant, Kim, & Kamionkowski 
(2004, hereafter EKK04) very recently showed that the total power generated by DF within 
the cooling radius is comparable to radiative loss from that region if the average mass-to-light 
ratio of galaxies is about 10. They also noted that this mechanism is self-regulating in the 
sense that gas would not be overheated since the DF of galaxies becomes ineffective when it 
attains high enough temperature for which galaxy motions become subsonic (see also Miller 
1986). 

While the results of the aforementioned work suggest that the total power supplied by 
the motions of member galaxies is comparable to X-ray luminosity of a typical rich cluster, 
it is sill uncertain whether it can balance the radiative cooling locally as well. If the DF of 
galaxy motions is to serve as a primary heat supplier in clusters in equilibrium, observed 
density and temperature profiles of clusters should be a natural consequence of the local 
heat balance between radiative cooling and DF-induced heating. Since optically-thin, X- 
ray emitting gas is prone to thermal instability (e.g., Fabian 1994), it is also an interesting 
question whether the DF-induced heating is able to suppress thermal instability completely 
or at least lengthen its growth time to the level comparable to the ages of clusters; even 
an equilibrium cluster would otherwise be subject to a mass dropout at its center (Kim & 
Narayan 2003a). 

In this paper, we take one step further from Miller (1986) and EKK04 to investigate 
equilibrium cluster models in which DF-induced heating is a main heating mechanism. We 
will first construct density and temperature distributions of the gas in clusters that are in 
strict hydrostatic and thermal equilibrium and compare them with observed profiles. We note 
that recent X-ray data from BeppoSAX, Chandra, and XMM-Newton observations indicate 
that gas temperature in rich clusters rapidly increases with radius for r £ (0.05 — 0.08)ri 80 , 
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remains roughly isothermal up to r ~ 0.2ri 8 o, and then gradually declines farther out (Ettori 
et al. 2000; De Grandi k Molendi 2002; Piffaretti et al. 2005; Vikhlinin et al. 2005, see also 
Markevitch et al. 1998 for ASCA results). Here, t is a virial radius the interior of which has 
the mean density of gas equal to 180 times the critical density. While declining temperature 
profiles at large radii may provide important clues to cosmological structure formation, they 
appear not to be directly relevant to the cooling flow problem. This is because gas density 
in clusters decreases as or faster than ~ r~ 2 at the outer parts, so that the radiative cooling 
time in the r ^ 0.2r 180 regions are much longer than the Hubble time. Indeed, the results of 
numerical simulations for cluster formation show that the shapes of declining temperature 
profiles are essentially independent of the presence of radiative cooling and/or supernova 
feedback (e.g., Loken et al. 2002; Motl et al. 2004). For this reason, when we compare our 
model results with observations, we pay attention to density and temperature profiles only 
in cooling regions with r <, 0.2r 180 (typically ~ 0.5 — 1 Mpc) beyond which the effect of 
radiative cooling is not serious and thus heating mechanisms are not required. 

The remainder of this paper is organized as follows: In §2 we evaluate the total heating 
rate resulting from the DF of galaxies using the formula given by Ostriker (1999) for gravita- 
tional drag force in a gaseous medium. While Miller (1986) and EKK04 approximated galaxy 
motions as being highly supersonic, we allow for both subsonic and supersonic motions of 
galaxies. In §3 we calculate equilibrium density and temperature profiles of galaxy clusters 
by assuming that the DF-induced heating is deposited at the locations of galaxies. Effects of 
distributed heating of DF and the radial mass variation of galaxies are discussed in §4. In §5 
we analyze local thermal stability of the gas while in §6 we present the results of numerical 
hydrodynamic simulations that investigate the evolution of initially out of equilibrium con- 
figurations. Finally, in §7 we conclude this work and discuss other potential consequences of 
DF of galaxies in clusters. 



2. Total Heating Rate 

Galaxies orbiting around the center of a galaxy cluster gravitationally induce wakes in 
the ICM. The gravitational interaction between the galaxies and wakes causes the former to 
lose their orbital kinetic energy and converts it into thermal energy of the ICM via either 
compressional heating (including shock) or turbulent dissipation of the gas kinetic energy 
in the wakes. Although the notion that DF of galaxies can generate heat in the hot ICM 
has been well recognized and widely studied, most of previous work on this subject assumed 
galaxies in highly supersonic motions, estimated only total heating rate, and compared it 
with total radiative loss rate in clusters (Ruderman & Spiegel 1971; Rephaeli & Salpeter 
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1980; Miller 1986; EKK04). However, the motions of galaxies are only slightly supersonic, 
with an average Mach number of about 1.5 (e.g., Sarazin 1988; Faltenbacher et al. 2004), 
and even become subsonic at the outer parts of clusters. In this paper we consider both 
subsonic and supersonic galaxy motions and adopt the general formula of Ostriker (1999) 
for the DF force in a gaseous medium. 

Following Ostriker (1999), we consider a galaxy of mass M g moving at velocity v through 
a uniform medium with density p and sound speed c s . The dynamical- friction force that the 
galaxy experiences is given by 



DF 



Anp(GM g fI 



v 3 



(1) 



where the efficiency factor / is defined by 



- \ |ln(l-^- 2 ) + ln(^/r min ), M > 1, V> 

with M. = v/c s being the Mach number of the galaxy motion and r min the effective size 
of a galaxy (Ostriker 1999). Note that for M. 3> 1, equations (1) and (2) with vt = r max 
corresponding to the maximum system size, become identical to the Chandrasekhar (1943) 
formula for the drag force due to collisionless particles. Although equation (1) is valid for 
a perturber moving on rectilinear trajectory through a uniform- density medium, numerical 
simulations show that it is a good approximation even in a radially-stratified, spherical 
system (Sanchez-Salcedo & Brandenburg 2001). 

We assume that the orbital velocities of N g galaxies that contribute to heating of the 
gas are described by the Maxwellian distribution, 

/(«) - j^fm^'^. <3) 

with a one-dimensional velocity dispersion oy. The total heating rate due to the DF is then 
given by 

( ^ =W .<_ F „. v) = <=£^/'\ ,4) 



c s \M, 

where the angular brackets denote an average over equation (3). 

Using equations (2) and (3), one can evaluate (I / M) numerically. Figure 1 plots {I / M) 
as functions of m = a r /c s for some values of In (vt/r niin ). For supersonic cases, the density 
perturbations in the wake are highly asymmetric and the regions influenced by the perturber 
shrink as m increases. This results in a smaller heating rate for larger m. In this case, one 
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can show (I/M) « (2/-n) 1 / 2 m~ 1 exp(— 0.5/m 2 ) In (vt/r min ) for m 3> 1. For a subsonically 
moving perturber, on the other hand, the perturbed density in the front and back sides of the 
perturber becomes more or less symmetric, producing a smaller heating rate with decreasing 
m. In the limit of vanishingly small m, (I/M.) — > m 2 . In general, (I/M.) peaks at m — 1, 
as Figure 1 indicates. 

Taking v t — 1 Mpc from a typical size of clusters and r min = 10 kpc as the galaxy size, 
In (vt/r min ) 4.6. Since m usually varies from 0.8 to 3, (I/M) w 2. Therefore, we have the 
total heating rate 

<*> - * «- - ® d^k) 2 (o^) t^r m • 

where n e and T denote the electron number density and gas temperature, respectively, and 
we adopt the solar abundance. Notice that the rate of gas heating given in equation (5) is 
similar to typical X-ray luminosity of rich clusters (e.g., Sarazin 1988; Rosati, Borgani, & 
Norman 2002), implying that there is a sufficient amount of energy available in the orbital 
motions of galaxies. This is essentially what led EKK04 to conclude that the dynamical- 
friction coupling between cluster galaxies and gas can provide thermal energy enough to 
compensate for radiative loss. However, it is still questionable whether DF of galaxies heats 
the gas in a right manner. That is, are the observed density and temperature profiles of the 
ICM a direct consequence of energy balance between heating by DF and radiative cooling? 
Is the intracluster gas heated by DF thermally stable? In what follows, we shall study 
models of clusters with DF-induced heating in detail by assuming hydrostatic equilibrium 
and thermal energy balance. 



3. Model With Localized Heating 

DF of galaxies is mediated by gravity which is a long-range force, so that heating 
of gas due to a single galaxy is likely to be well distributed throughout its wake. The 
functional form of heat distribution is unknown and its derivation may require numerical 
simulations, which are beyond the scope of the present work. In this section, we instead 
make a simplifying assumption that heat is deposited locally at the galaxy position. The 
effects of heat distribution will be discussed in §4. 
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3.1. Local Heating Function 

Consider a spherically symmetric distribution of galaxies, with the number density given 
by n g (r) at radius r. The local heating rate per unit volume due to DF is given by 

«r)> = AMGM ' ?{IIM) n,(r), (6) 

Cs 

where we assume equal galaxy mass. To find n g (r), we assume that the orbits of galaxies are 
isotropic and isothermal. The usual Jeans equation (e.g., eq. [4-55] of Binney & Tremaine 
1987) then reads 

where $ is the gravitational potential. 

Under the NFW distribution of dark matter 

_ M /2n 

PNFW " r(r + r s r (8) 

with a scale radius r s and a characteristic mass M , the gravitational acceleration is given 
by 



d$ 2GM 
dr r 2 s 



\n(l + x) 



, . ■ (9) 

where x = r/r s is the dimensionless radius (Navarro, Frenk, & White 1997; Klypin et al. 
2001). Combining equations (7) and (9) together, one can easily find 



gO*Q 
f g{x')d i x''' 



n g (x) = n a (0)g(x) = N g /JJ,^, , (10) 



where 

with the dimensionless parameter i] defined by 



(11) 



2GM 

7 i = ^r-^- (i 2 ) 



Note that r] measures the ratio of the gravitational potential energy of a galaxy in a cluster 
to that galaxy's kinetic energy at r = r s . 

While the functional form of g(x) looks strange, it behaves quite well near x ~ and 
gives density profiles similar to the isothermal /9-model of gas distributions (Makino et al. 
1998) or to the King model of the observed galaxy distributions (e.g., Girardi et al. 1998). 
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For example, Girardi et al. (1998) found that the best-fit distribution of galaxies in A1795 is 
n g{ r ) = n g(0)[l + ( r /-Rc) 2 ] -1 ' 27 with a core radius R c = 43 pc, which is plotted as a solid line 
in Figure 2. On the other hand, X-ray and optical observations indicate M = 6.6 x 1O 14 M , 
r s = 460 kpc, and a r « 800 km s" 1 for A1795 (Girardi et al. 1998; Ettori et al. 2002), 
corresponding to rj m 19. Figure 2 plots as a dotted line the g(x) curve with rj = 19, which 
is in good agreement with the observed galaxy distribution. 

We take 2r s as the outer boundary of a cluster. Inserting equation (10) into equation 
(6), one obtains 

4np(GM 9 fN 9 I I \ 
c s r s 3 / g{x')d 3 x> \A</ 

1.7 x lO^V-^^sooM^^eT^^//^)^^) erg cm" 3 s" 1 , (13) 

where r M6 o = r s /(460 kpc), N gt5O0 = N g /500, M fl)11 = M g /(10 U M G ), and T kcV is the 
temperature of the gas in units of keV. Equation (13) is our desired equation for the volume 
heating rate due to the DF of galaxies. Note that (I / M) ~ 2 at m ~ 1 and depends on 
temperature. 



(e(x)) = 



3.2. Equilibrium Model 

We look for cluster density and temperature distributions in which the hot gas satisfies 
both hydrostatic equilibrium and energy balance between radiative cooling and DF-induced 
heating. For thermal bremsstrahlung, the rate of energy loss per unit volume is given by 

j = 7.2 x 10- 2A n 2 e T^ erg cm" 3 s _1 , (14) 

(Rybicki & Lightman 1979; Zakamska & Narayan 2003). From equations (13) and (14), the 
condition of thermal energy balance yields 

n e T keV = 2.4 x 10- 2 r-l m N g>50Q Ml 11 (I/M)g(x). (15) 

Neglecting the weak temperature dependence of (I/A4) and assuming that M g does not 
vary with radius, equation (15) states that the gas pressure in equilibrium should trace the 
distribution of galaxies. 

Hydrostatic equilibrium of the gas requires 

c\ d(n e T kcV ) _ d<& ^ 
n e dr dr 1 
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where Co = 390 km s 1 is an isothermal sound speed at Tk e v = 1- Substituting equations (9) 
and (15) into equation (16) and using equation (11), we obtain 



■2 



^-\n(I/M) = - V (-£--1 



ln(l + x) 



x 2 x(l + x) 



(17) 



Since (I/A4) depends on T ke y, equation (17) can be integrated as an initial value problem 
to yield T keV (x). The equilibrium profile for n e (x) can then be found from equation (15). 
Since n e (x) depends rather sensitively on less- well known quantities such and N g and M 2 
and is easily influenced by the mass profile of galaxies, we concentrate on T ke y(x) that is 
independent of the properties of galaxies other than ay. 

Equation (17) has a trivial isothermal solution 

T >» = 65 (lo^) 2 keV - < 18 > 

independent of the radius. The corresponding density profile is 

n e = 7.4 x 10- 3 r- 4 3 60 a- 3 2 iV Si500 M 9 2 11 ( 7 (x), (19) 

where a rj3 = oy/(10 3 km s _1 ), indicating that the gas density follows the galaxy distribution 
exactly. Considering uncertainties in the values of r Sj460 , o^ioooj N 9t500 , and M^n, 1 equation 
(19) gives n e ~ 0.01 — 0.2 cm~ 3 at the cluster centers, which is not much different from 
observed values. Note that the Mach number corresponding to Ti SO is m = r )~ l l 2 . For later 
purposes, we define the transonic temperature 

-^sonic = -^isoj (20) 

7 

corresponding to unity Mach number of galaxy motions. 

The isothermal solution (18) is consistent with what has long been known as the scaling 
relation between the gas temperature and the velocity dispersion of galaxies for many clusters 
(e.g., Mushotzky 1984; Edge & Stewart 1991; Lubin k Bahcall 1993; Wu et al. 1999, see 
also Sarazin 1988; Rosati, Borgani, & Norman 2002 for review). It explains why the values 
of the (3 (= nm p a 2 /kT) parameter in the f3 models for gas distributions are close to unity 
(Lubin & Bahcall 1993; Bahcall & Lubin 1994; Girardi et al. 1996b, 1998). Physically, the 
scaling relation implies that the plasma and the galaxies are well relaxed under the common 
gravitational potential. It is uncertain whether the observed scaling law results truly from 



1 For example, the average mass of galaxies would increase up to ~ 5 x 10 11 M Q if the contribution from 
dark halos is taken into account (e.g., Zentner & Bullock 2003). 
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the DF coupling between the ICM and galaxies, but equation (18) suggests that the latter 
certainly makes the former tighter. 

Although gas in the regions surrounding cooling cores of clusters has nearly constant 
temperatures close to the virial values, recent X-ray observations exhibit positive radial 
gradients of gas temperatures at the central ~ 300 kpc regions (Peterson et al. 2001, 2003), 
while showing declining temperature distributions at far outer regions (Markevitch et al. 
1998; De Grandi & Molendi 2002; Piffaretti et al. 2005; Vikhlinin et al. 2005). Thus, the 
isothermal solution cannot describe the observed temperature distributions for the whole 
range of radii. We are particularly interested in the regions within and adjacent to the 
cooling radius in which gas density is high enough to experience significant radiative cooling. 
To check whether the general solutions of equation (17) produce temperature and density 
distributions similar to the observed in such regions, we solve it numerically. We adopt 
v = 19, r s = 460 kpc, N g = 500, M g = 10 n M o , and a r = 10 3 km s" 1 corresponding to 
T iso = 6.5 keV. We choose a value for the central temperature T(0) and then integrate 
equation (17) from r = to 2r s . The resulting temperature and density profiles for a few 
selected values of T(0) are plotted in Figure 3. 

The spatial behavior of the solutions depends critically on T(0). When T(0) > T sonic , 
galaxy motions are subsonic and d(I / M) / dm > from Figure 1. If T(0) > T; so (or m < 
7~ 1/2 ), equation (17) gives d{I/M)/dr = -(m/2T)(d(I/M)/ dm) (dT/dr) > for T(0) > 
T iso , implying dT /dr < 0. As temperature monotonically decreases with radius, m increases 
until it reaches the isothermal value 7~ 1 / 2 where T = T iso . One can similarly show that 
dT/dr > for T iso > T(0) > T sonic , and thus temperature slowly increases toward T iso . As 
long as T(0) > T sonic , the solutions asymptote to the isothermal ones (eqs. [18] and [19]) as 
the radius increases. On the other hand, dT/dr < when T(0) < T sonic (or m > 1). Since 
m increases with decreasing T and since d{I/M)/dr does not change its sign for m > 1, 
dT/dr < is satisfied for all radii. The decrease of temperature is much faster than that 
of g(x), resulting in a unrealistic distribution of electron number density that increases with 
radius. 

As Figure 3 shows, the local heat balance with DF-induced heating yields rising temper- 
ature profiles, if and only if T sonic < T < T iso . While this is a tantalizing result considering 
the central depression of temperatures seen in cooling-flow clusters, the required range of 
temperatures is very narrow. If we identify T; so with the virial temperature, the central 
temperatures in our equilibrium models must be larger than T sonic = 0.6T iso (for 7 = 5/3), 
which is about twice larger than the observed values of T(0) ~ (0.3 — 0.4)T iso (e.g., Peterson 
et al. 2001, 2003; Allen et al. 2001). By having too stringent upper and lower limits of 
temperatures, therefore, heating by DF of galaxies alone is unlikely to explain observed tern- 



perature and density distributions of galaxy clusters. We note however that the tight range 
of temperatures for the rising temperature profiles may be caused by the basic assumptions 
of this section, namely that DF-induced heating of the ICM is all localized at the positions 
of galaxies, and that all galaxies have equal mass. We relax these assumptions in the next 
section. 



4. Effects of Distributed Heating 

In the presence of viscosity and/or turbulent magnetic fields in the ICM, the kinetic 
energy lost by a galaxy via DF will eventually be converted into heat, rather than bulk 
motions distributed onto the gravitational wakes of the galaxy. Deiss & Just (1996) used a 
quasi-linear fluctuation theory to derive the spatial structure (in Fourier space) of mechanical 
heating in a turbulent medium self-consistently driven by DF of many galaxies. To our 
knowledge, there is no published work that study spatial heat distribution (in real space) 
caused by DF. If ICM heating by DF of galaxies occurs through turbulent dissipation of 
gas kinetic energy, finding the functional form of the heat distribution would not be viable 
unless the characteristics of ICM turbulence and related processes are prescribed (Deiss, 
Just, & Kegel 1990). Instead of attempting to derive a realistic heat distribution, in this 
work we simply adopt a Gaussian function to parametrize the spatial extent to which heat 
is distributed 2 . Our aim in this section is to examine the effects of heat distribution on 
equilibrium structures in comparison with the localized heating models. Since masses of 
galaxies inside a cluster are likely to vary with the radius within the cluster, we also allow 
for spatially varying galaxy masses. 

Assuming that heat distribution follows a Gaussian profile centered at the location of 
the galaxy, we write the heating rate per unit volume as 

(e(r)) = h(r), 

c s 

where 

h(r) = ^fj n g {r')M g {r'fe-^' l ldr\ (22) 

is a convolution of n g M^ with the Gaussian function with scale length lh- In writing equations 
(21) and (22), we assume that density and temperature of the gas do not change significantly 



2 Wc also tried with logarithmic heat distribution functions used in EKK04 and found that results are 
similar to those with Gaussian functions presented in this section. 
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inside the wake of a galaxy, which is acceptable within the framework of equation (1). The 
condition of hydrostatic equilibrium and thermal energy balance is then reduced to 



d 

dx 



\n(I/M) = 




ln(l + x) 

x 2 



x(l + x) 



1 



dha.h(x) 



dx 



(23) 



It can be easily verified that equation (23) recovers equation (17) in the limit of 4 — > 0. 

Figure 4 plots the solutions of equation (23) for the case of constant M g . More widely 
distributed heating is equivalent to localized heating with a flatter distribution of galaxy 
number density, causing the \d\nh/dx\ term in equation (23) to be smaller. Consequently, 
an equilibrium cluster with a larger value of lh tends to make (I/Ai) decreases faster with 
increasing r. This implies that for T(0) < T sonic (T(0) > T sonic ), the temperature in models 
with non-zero lh decreases (increases) more rapidly than in the lh — counterpart, although 
the difference between models with lh/r s = 0.1 and 1.0 is negligible for small r. Note 
that even when T(0) > T iso , for which localized-heating models have declining temperature 
profiles, equilibrium temperature in distributed-heating models increases in the regions with 
r < l h and eventually converges to T iso at r £ (2 — 3)^. We see that distributed heating 
does not change the condition T son i C < T < Tj SO for the existence of rising temperature 
distributions. The distributed-heating models do not work better than the localized-heating 
model in terms of producing temperature distributions similar to observations; it in fact 
aggravates the situation by making all of the gas in clusters nearly isothermal except only 
in the central ~ 1 kpc regions. 

Another factor that may affect the equilibrium structure is the radial mass variation of 
cluster galaxies. Although galaxies located near the central parts of clusters tend to have 
a large fraction of luminous matter, (e.g., Biviano et al. 2002; Mercurio et al. 2003), the 
strong tidal stripping of dark-matter halos is likely to cause the total (luminous + dark) 
masses of individual galaxies to be smaller toward the cluster center (Zentner & Bullock 
2003; EKK04; Nagai & Kravtsov 2005). Motivated by this consideration, we adopt several 
algebraic forms of M g (r) that either monotonically decrease or increase by about a factor of 
20 or less from the cluster center to r = 2r s , and then solve equation (23). The resulting 
temperature profiles are almost identical to those of the constant-mass cases presented in 
Figure 4. Because the distribution of galaxy masses appears logarithmically in equation 
(23), it does not considerably affect temperature structures, although it dramatically changes 
electron number-density profiles. 
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5. Thermal Stability 



We now discuss thermal stability of a system in which radiative cooling is locally bal- 
anced by DF-mediated heating. For spatially localized heating, the net loss function pC is 
given by 

pC = j — (e) = ap'T 1 ' 2 - PpT-W+r, (24) 

where the power index p accounts for the local temperature dependence of (I/A4) in a 
piecewise manner, and the positive constants a and (3 contain all information on the cluster 
properties other than gas density and temperature. Figure 5 plots p = d\n(I / M) / d\nT 
as a function of m, which gives p ~ at around the transonic temperature (m fa 1), while 
p fa 0.3 — 0.5 for supersonic temperatures (m > 1). The curves asymptote to —1 for m <C 1 
and to 0.5 for m ^> 1, as the asymptotic formulae given in §2 suggest. Note that p is 
insensitive to the choice of \n(vt/r min ) as long as m > 0.8. 

Thermal stability of a system can be easily checked by using the generalized Field 
criterion which states that the system is thermally unstable to isobaric perturbations if 



d{C/T) 


1 


'd(pC) 


dT 


P pT 


dT 



P d{pC) 
T dp 



-p 



ap 

f3/2 



<0, 



(25) 



where the second equality assumes pC = corresponding to thermal equilibrium (Field 1965; 
Balbus 1986). Since p > for m > 1 (and p < for m < 1), this implies that the ICM 
is thermally unstable (stable) if heating is caused preferentially by galaxies moving at a 
supersonic (subsonic) speed. Dense inner parts of galaxy clusters, where the cooling time is 
less than the Hubble time, are filled with low-temperature gas such that galaxy motions are 
readily supersonic, implying that DF-induced heating is unable to quench thermal instability 
in those cooling cores. Whether thermal instability has important dynamical consequences 
on cluster evolution can be judged by considering its growth time which amounts to 
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(26) 



(cf. Kim & Narayan 2003a). Even if cooling cores are thermally unstable, therefore, the 
virulence of thermal instability may or may not manifest itself during the lifetime of clusters 
depending on the value of p. For example, for a cluster with a r = 10 3 km s" 1 and T = 2 
keV at the central regions, m = 1.4 and p = 0.23 from Figure 5, giving t grow — 4.2 Gyr. 
This is almost comparable to the ages of clusters since the last major merger event (~ 7 Gyr 
for massive clusters; Kitayama & Suto 1996), suggesting that thermally instability may be 
dynamically unimportant for practical purposes. 
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The inability of DF-induced heating to suppress thermal instability seems not to be a 
fatal problem as long as the following two conditions are met: i) cooling cores are in thermal 
equilibrium, a basic premise of the stability analysis; and ii) galaxies near the central parts 
move at near-transonic speeds (m £2), so as to have a sufficiently small value of p. While 
the second condition appears to be easily satisfied in clusters, the results of the preceding 
section suggests that the first condition may not be so if DF-induced heating is to balance 
radiative cooling for the whole range of radii. Since the cooling time outside the cooling 
radius is longer than the Hubble time, one may argue that material beyond the cooling 
radius does not need any heating source and that it is sufficient for DF of galaxies to supply 
heat to gas only in cooling cores. We address this possibility in the next section. 

As a final remark of this section, we compare our results with previous work. Formally, 
equation (25) implies that gas subject to DF-induced heating has a critical temperature 
above (below) which it becomes thermally stable (unstable) and that this critical temper- 
ature corresponds exactly to unity Mach number of galaxy motions. This result is quite 
different from those of previous work which made some approximations on the heating func- 
tion. For instance, Miller (1986) argued that gas with drag heating is always thermally 
unstable and rapidly turns into a multi-phase medium, a consequence of his neglecting the 
temperature dependence of the heating function (corresponding to p = 1/2 in eq. [24]). 
Although Bregman & David (1989) found a similar critical gas temperature, they consid- 
ered only the supersonic regime and used an approximate heating function from the shock 
jump conditions, which makes their critical temperature smaller than ours by a factor of 
about three. Our result, based on the general formula of Ostriker (1999) for the drag force 
applicable for both supersonic and subsonic cases, represents the situation better. 



6. Time Dependent Approach 

We have shown in the preceding sections that an equilibrium cluster subject to DF- 
mediated heating has a temperature profile that either decreases (for T > T iso or T < T sonic ) 
or increases (for T sonic < T < T iso ) with radius. While the radially increasing temperature 
profile is attractive, the requirement that thermal energy is balanced locally in the entire 
regions out to 2r s causes the central gas temperature to be no less than 0.6 times the 
virial temperature, about a factor two larger than observations. Since radiative cooling is 
important only in the dense central parts inside the cooling radius, however, it would be 
sufficient for DF of galaxies to provide heat only in such a cooling core instead of the entire 
regions. In addition, thermal instability would not be an issue provided that gas is in thermal 
equilibrium and that galaxy motions are near transonic (§5). Therefore, it may be possible 
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for a cluster to start from an arbitrary non-equilibrium state and evolve slowly under both 
radiative cooling and DF-mediated heating, ending up with an equilibrium core in which 
thermal instability has a very long growth time. The resulting temperature profile, albeit 
consistent with observations, need not represent an equilibrium for the whole range of radii. 

To test this idea, we have run a number of one-dimensional hydrodynamic simulations 
in spherical polar coordinates. We set up a logarithmically spaced radial grid, with 400 
zones, from 1 kpc to 1 Mpc, and explicitly implement the net cooling-heating function 
described by equation (24). We also include the effect of gravitational potential due to 
passive dark matter using equation (9) with Mq = 6.6 x 10 14 M Q and r s = 460 kpc; the 
DF-coupling between gas and the dark matter and back reaction on the latter component 
are not taken into account. As initial conditions, we consider spherically symmetric clusters 
and take n e = n e (0)g(x) cm -3 with the central electron density n e (0) = 0.002, 0.005, 0.01, 
0.02 cm -3 , and fix T = 6.5 keV independent of the radius for all models. We adopt r s = 460 
kpc, N g = 500, and M g = lO n M , and consider both localized heating and distributed 
heating with 4/r s =0.1 and 1. Note that the corresponding equilibrium density profile has 
n e (0) = 0.0074 cm" 3 when heating is localized. We adopt the outflow boundary conditions 
for scalar variables (i.e., density, energy, etc.) that assign the same values in the ghost zones 
as in the corresponding active zones. For the radial velocity, we allow it to vary as a linear 
function of radius at the inner boundary, while fixing it to be zero at the outer boundary 
(e.g., Kim & Narayan 2003a). Using the ZEUS hydrodynamic code (Stone & Norman 1992), 
we solve the basic hydrodynamic equations and follow the nonlinear evolution of each model 
cluster. The resulting radial profiles of electron number density and temperature of model 
A with n e (0) = 0.02 cm" 3 and model B with n e (0) = 0.005 cm -3 , both of which assume 
spatially localized heating, are shown in Figure 6 for a few time epochs. 

Model A, which initially has larger density than the equilibrium value everywhere, im- 
mediately develops radial mass inflows. As time evolves, the temperature drops and the 
density increases. Compared to pure cooling cases which fairly well maintain near isobaric 
conditions (e.g., Kim & Narayan 2003a), DF-induced heating is found to cause thermal pres- 
sure to increase with density as P oc p a3 . Since j oc T~ 2A and (e) oc T~ L9+P with p > for 
supersonic galaxy motions, however, the system is thermally unstable and cooling occurs at 
a much faster rate as temperature decreases. Model A experiences a cooling catastrophe in 
less than 4 Gyr. On the other hand, model B with initial overheating becomes hotter and 
less dense, which in turn tends to increase the ratio of heating rate to the cooling rate for 
p ?t — 0.5. As the gas temperature increases, the Mach number of galaxy motions become 
smaller, reducing the value of p. The cooling rate will therefore eventually exceed the heating 
rate and the cluster evolution will be reversed. In the case of model B, this turnaround will 
occur at t — 70 Gyr when T ~ 18 keV (or m ~ 0.5) is reached. 
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The evolution of models with different density and different length scales for heat dis- 
tribution is qualitatively similar to those of models A and B, namely that clusters catas- 
trophically evolve toward vanishing central temperature if cooling dominates initially, while 
clusters with initial excessive heating heat up steadily. This suggests that DF-induced heat- 
ing does not naturally lead non-equilibrium clusters to thermally stable, equilibrium cores. 
As shown in §3 and §4, DF of galaxies does not explain the observed temperature distribu- 
tions of clusters if the condition of thermal energy balance is imposed for all radii. Although 
it is enough for DF-induced heating to balance radiative cooling only in cooling cores, such 
cooling cores do not naturally form from non-equilibrium states. Unless the properties and 
galaxies and ICM are fine tuned, small departures from an equilibrium state rapidly evolve 
into an extreme configuration. Therefore, we conclude that DF-induced heating alone is not 
likely to account for the absence of cold gas in the centers of galaxy clusters. 

Although heating by DF from galaxies does not appear to provide a complete solution to 
the cooling flow problem, we see that the DF-induced heating can still offset a considerable 
amount of radiative cooling. The isobaric cooling time is given by 

7 ( P \ _ 0-96 / n e \-i / kT 

COQl -i-l\pC) l-ft/(n e T keV ) V0.05cm- 3 7 V 2 keV 

where 1Z = 0.024r~4 60 A r 9i500 M| 11 (//.M) represents the contribution of the DF-induced heat- 
ing. Figure 7 plots as solid lines the cooling time as a function of n e for N gj500 M^ u = 1 or as 
a function of N g ^M gll with varying density. The cases without any heating are compared 
as dotted lines. The time epochs when clusters experience cooling catastrophe in the numer- 
ical simulations are marked in Figure 7a as solid circle, triangle, cross, and open circle for 
no heating, the distributed heating with l^/r = 1.0 and 0.1, and the localized heating cases, 
respectively. The numerical results are in good agreement with the theoretical prediction. 
Models with distributed heating tend to have longer cooling time. When N gi500 M gll = 1 
and n e (0) = 0.05 cm -3 , Figure 7a shows that DF-induced heating lengthens the cooling 
time by about 13% compared to the pure cooling case. However, as Figure 7b indicates, the 
offset of cooling by DF-induced heating is increasingly larger as N gi500 M gll increases. The 
increment of the cooling time could be as large as 140% when n e (0) ~ 0.03 — 0.1 cm -3 and 
Ng£ooMg t n ~ 3 — 7, suggesting that thermal energy supplied by DF of galaxies is by no 
means non-negligible. 

7. Conclusions and Discussion 

Friction of galaxy motions via the gravitational interaction with their own gravitation- 
ally induced wake in the ICM has often been invoked as an efficient heating mechanism of 
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the gas (e.g., Miller 1986; Just et al. 1990; Deiss & Just 1996; EKK04). This idea of DF- 
mediated heating of the ICM is quite attractive because there is sufficient energy available 
in galaxy motions and the mechanism is self-regulating; it operates effectively only when the 
temperature of gas is in a certain range, which happens to be the typical gas temperatures 
in the cooling cores of galaxy clusters. In this paper, we take one step further from Miller 
(1986) and EKK04 to calculate equilibrium density and temperature profiles of the hot gas 
heated by DF of galaxies. Instead of restricting to cases where galaxy motions are all su- 
personic, we use the general formula derived by Ostriker (1999) for the drag force that takes 
allowance for both supersonic and subsonic galaxy motions in a gaseous medium. We show 
that the total heating rate due to the DF of galaxies in a typical rich cluster is comparable 
to its total X-ray luminosity, confirming the results of Miller (1986) and EKK04 (see §2). 

Next, we derive the local heating function (eq. [13]) assuming that the orbits of galaxies 
are isotropic and isothermal under the NFW distribution of dark matter and that the kinetic 
energy lost by a galaxy is deposited into heat at the location of the galaxy. The condition that 
the gas is in hydrostatic equilibrium and maintains energy balance requires the temperature 
profile to be one of the following three kinds: isothermal, with T = Tj SO ; a decreasing 
profile with radius when T < T sonic or T > T iso ; an increasing profile when T sonic < T < T[ so , 
where T iso and T sonic denote the temperatures corresponding to unity Mach number of galaxy 
motions with respect to the isothermal and adiabatic sound speeds, respectively (eqs. [18] 
and [20]). The isothermal solution is interesting because it quite well describes the observed 
scaling relationship among clusters properties. Although the rising temperature solution is 
attractive since clusters usually show a temperatures drop at the central regions, it strictly 
requires the central temperature to be no lower than 0.6 times the virial temperature, which 
is roughly twice smaller than the observed values (§3). We also consider cases in which 
DF-induced heating is distributed in space according to a Gaussian form, and/or the masses 
of cluster galaxies vary over radius, and find that the stringent limit of T > T sonic = 0.6T iso 
for rising temperature distributions remains unchanged (§4). 

Using the local heating function we have derived, we analyze thermal stability of the 
gas subject to radiative cooling and DF-mediated heating (§5). When galaxy motions are 
subsonic, the heating rate that deceases steeply with temperature suppresses thermal insta- 
bility completely. On the other hand, supersonic galaxy motions for which the heating rate 
is relatively insensitive to temperature are unable to stop the growth of thermal instability. 
The growth time in the presence of DF-induced heating is at least twice that in the pure 
cooling case, and becomes progressively longer as the Mach number of galaxies decreases. 
When galaxies move at slightly supersonic velocities with Mach number less than 2, which 
is very likely at cluster centers, the growth time becomes comparable to the ages of clusters. 
This implies that thermal instability, even if operates, may have insignificant dynamical 
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consequences on cluster evolution. 

Noting that regions outside the cooling radius need not be in strict thermal equilibrium, 
we look for a possibility using numerical hydrodynamic simulations that clusters evolve 
from an arbitrary non-equilibrium state and form current cooling cores in which DF-induced 
heating balances radiative cooling (§6). We find that clusters that were initially dominated 
by cooling unavoidably develop radial mass inflows and decreases their central temperatures 
in a runaway fashion, whereas clusters with initial overheating slowly heat up and result 
in radially decreasing temperature profiles. Equilibrium solutions therefore do not appear 
to form an attracting set for galaxy and gas configurations, suggesting that when DF from 
galaxies is the sole heating source, it is extremely difficult to obtain equilibrium cores by 
smoothly evolving non-equilibrium clusters (even if in some cases the cooling catastrophe 
can be deferred so as to occur on a longer timescale). 

Putting together all the results of this paper we conclude that DF of galaxies alone, 
albeit an interesting heating mechanism with a lot of available energy, cannot be the lone 
heating agency to balance radiative cooling in rich galaxy clusters. We nonetheless note that 
the heating due DF of galaxies could considerably lengthen the cooling time, depending on 
the value of N g Mg, and thus should not be neglected in energetics of galaxy clusters. 

One of the key assumptions made in this paper is that all the kinetic energy lost by 
galaxies via DF is transferred to the thermal energy of the surrounding gas. The rationale 
behind this assumption is that the superposition of the wakes produces turbulence in the 
ICM and the kinetic energy injected into the ICM turbulence at saturation cascades down 
along the Kolmogorov-like energy spectrum, and eventually transforms into heat through 
viscous dissipation at small scales. Another possibility is that a large fraction of the kinetic 
energy of galaxies is used to merely enhance the level of the turbulence instead of being 
converted into heat. This may occur when the turbulence is not fully developed yet. 

The characteristics of the ICM turbulence is not well known, yet observations and 
numerical simulations suggest an average velocity dispersion i>t U rb ~ 200 — 400 km s" 1 on 
scales of A ~ 5 — 20 kpc (Ricker & Sarazin 2001; Carilli & Taylor 2002; Churazov et al. 2004; 
Schuecker et al. 2004; Faltenbacher et al. 2004). The associated turbulent kinetic energy 
is (l/2)M gas t; t 2 urb ~ 9 x 10 61 erg for the total mass M gas ~ 1O 14 M of the ICM in a rich 
cluster. Let us assume that this energy is supplied solely by DF of galaxies at a rate given 
by equation (5). If the level of turbulence keeps increasing without dissipation, it would take 
about 7 Gyr for the DF of galaxies to feed the ICM turbulence to the observed level, which 
is almost comparable to the lifetime of the cluster. Conversely, if the turbulence is fully 
developed and in a steady state such that the energy injection rate by DF is equal to the 
rate of dissipation, it would have iWb — {2XE/M g 

as) ^ ~ 60 km s * for A ~ 20 kpc, which is 
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too small to explain the observations. All of these imply that the contribution from the DF 
of galaxies to the ICM turbulence is not considerable (see also Sanchez-Salcedo et al. 1999). 
If the dissipation of the (fully developed) turbulence is to provide enough heat to balance 
radiative cooling, as suggested by Churazov et al. (2004) and Fujita et al. (2004), the energy 
injection should not be entirely due to the DF of galaxies; it requires other energy sources, 
including jets from active galactic nuclei and mergers with smaller groups or clusters. 

The assertion that the conversion of galaxy kinetic energy into thermal energy of the 
gas may not be so effective is supported by the results of Faltenbacher et al. (2004) who 
studied using numerical simulations cluster formation in ACDM cosmology, with the DF 
of galaxies included explicitly. They found that (1) the motions of galaxies at the present 
epoch are slightly supersonic, with an average Mach number of M. ~ 1.4; (2) gas within 
the virial radius has a velocity dispersion t> turb ~ (0.3 — 0.5)c s resulting probably from infall 
motions of galaxies and small groups, with small contribution from the DF of galaxies; 
and (3) the clusters still suffer from the cooling catastrophe. Although higher-resolution, 
more-controlled simulations are required to make a definitive statement, the last point of 
their results suggests that the system is thermally unstable or heating by DF may be quite 
inefficient — for, as we have seen in the previous section, systems evolved from arbitrary 
initial conditions do not generally tend to the equilibrium solutions, which do not therefore 
form an attracting set; and for some initial parameters the delay in the cooling catastrophe 
that the DF mechanism leads to is not significant. 

We finally note that several issues of potential importance were not investigated in this 
paper. One of these relates to the galaxy velocity anisotropy. It is easy to show that, for 
isothermal distribution with constant anisotropy, it is possible to obtain equilibrium solu- 
tions with smaller central gas temperatures, more in line with observations. However these 
tend to have unrealistic galaxy number density distributions. The equilibrium and stability 
of gas configurations in more realistic anisotropic models have not been investigated. We 
have adopted the general formula of Ostriker (1999) for the drag force on a galaxy moving at 
an arbitrary speed. This is an improvement over previous studies (e.g., Miller 1986; EKK04) 
that considered only supersonic galaxy motions. While its explicit dependence on the Mach 
number enabled us to explore temperature profiles of equilibrium clusters, the formula still 
assumes a galaxy moving in a straight-line orbit through a uniform gaseous medium. Clearly, 
the ICM is radially stratified and galaxies follow curved rather than rectilinear trajectories. 
In a collisionless system, Just & Peharrubia (2004) recently showed that the density inho- 
mogeneity reduces the Coulomb logarithm of the Chandrasekhar formula by limiting the the 
maximum impact parameter to the local scale length of density variation (see also Hashimoto 
et al. 2003; Spinnato et al. 2003), and induces an additional drag force in the lateral direc- 
tion of the galaxy motion (see also Binney 1977). They further showed that the reduction 
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of the Coulomb logarithm makes the orbital energy loss 15% less effective, inhibiting the 
circularization of the orbit, while the additional tangential drag force has a negligible effect 
on the orbital evolution. Similar effects are expected to occur in a gaseous medium, yet their 
consequences are not known quantitatively. A related important question is how heat (or 
turbulent) energy is distributed in the wakes that form in the turbulent, inhomogeneous, 
magnetized ICM. A real assessment of the effects of DF to dynamical evolution of galaxy 
clusters requires answers to the above questions, which will direct our future research. 
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Fig. 1. — Plots of (I / M.) as functions of m = a r /c s . Various curves corresponds to 
In (vt/r min ) = 4, 4.6, 6, 8, 10 from bottom to top. 
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Fig. 2. — Comparison of g(r) with i] = 19 (dotted) with the observed distribution (solid) of 
galaxies in A1795. The core radius is taken to be R c = 43 pc (Girardi et al. 1998). 
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Fig. 3. — (a) Equilibrium temperature and (b) electron density profiles of clusters in which 
DF-induced heating is spatially localized. The chosen parameters are rj = 19, r s = 460 kpc, 
N g = 500, M g = fO n M and a r = fO 3 km s -1 . Dotted and dashed lines in (a) respectively 
indicate the characteristic isothermal solution T; so = 6.5 keV and the transonic temperature 
T sonic = T iso /7 = 3.9 keV at which the Mach number of galaxy motions is unity. The density 
distribution corresponding to the isothermal solution is drawn as a dotted line in (b). Note 
that temperature increases with radius only when T sonic < T(0) < T iso . 
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Fig. 4. — Equilibrium (left) temperature and (right) density profiles of cluster models with 
distributed heating, which is modeled by a Gaussian function with scale length l h . The 
chosen parameters are the same as in Figure 3. When T(0) = 2, 4 keV, temperature profiles 
with lh/r s = 0.1 and 1.0 are identical to each other. See text for discussion. 
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Fig. 5. — Logarithmic slope p = d\n(I / A4) / dlnT = — 0.5d\n(I / A4) / dlnm as a function of 
Mach number m for In (vt/r min ) = 4, 4.6, 6, 8, 10 from top to bottom. The curves asymptote 
to — 1 for m 1 and to 0.5 for m > I. With positive values of p, DF-induced heating by 
supersonically moving galaxies with m > 1 is unable to suppress thermal instability. 
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Fig. 6. — Evolution of electron number density (top) and temperature (bottom) for (left) 
cluster model A with n e (0) = 0.02 cm~ 3 and (right) model B with n e (0) = 0.005 cm~ 3 . The 
initial temperature is 6.5 keV for both models. Model A, which initially has a larger density 
than the equilibrium value, inescapably experiences a cooling catastrophe, while model B 
with lower initial density steadily heats up. 
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Fig. 7. — Cooling time scale at the cluster center (a) as a function of n e for fixed N gt500 Mg u = 
1 and (6) as a function of N g M g for n e = 0.1,0.05,0.01 cm' 3 . Solid line is for the case 
when spatially localized, DF-induced heating is considered, while dotted line corresponds to 
the case without any heating. Various symbols in (a) represent simulation outcomes, with 
solid circle, triangle, cross, and open circle referring to no heating, distributed heating with 
h/fs — 1-0 and 0.1, and localized heating cases, respectively. 



